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ABSTRACT 

We present a X-ray spectral analysis of a large sample of 25 'bare' active galactic nuclei, 
sources with little or no complicating intrinsic absorption, observed with Suzaku. Our work 
focuses on studying the potential contribution from relativistic disc reflection, and examining 
the implications of this interpretation for the intrinsic spectral complexities frequently dis- 
played by AGN in the X-ray bandpass. During the analysis, we take the unique approach of 
attempting to simultaneously undertake a systematic analysis of the whole sample, as well as 
a detailed treatment of each individual source, and find that disc reflection has the required 
flexibility to successfully reproduce the broadband spectrum observed for all of the sources 
considered. Where possible, we use the reflected emission to place constraints on the black 
hole spin for this sample of sources. Our analysis suggests a general preference for rapidly 
rotating black holes, which if taken at face value is most consistent with the scenario in which 
SMBH growth is dominated by prolonged, ordered accretion. However, there may be obser- 
vational biases towards AGN with high spin in the compiled sample, limiting our ability to 
draw strong conclusions for the general population at this stage. Finally, contrary to popular 
belief, our analysis also impUes that the dichotomy between radio loud/radio quiet AGN is 
not solely related to black hole spin. 
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1 INTRODUCTION 

Knowledge of the spin distribution for the supermassive black 
holes powering active galactic nuclei (AGN) provides important 
information on the growth history of these black holes, which 
would otherwise be difficult to observationally constrain. The ac- 
cretion and merger events experienced by these nuclear black 
holes over their lifetimes can impart enough angular momentum 
to s ubstantially modify the black hole spins from t heir natal val- 
ues jModerski & Sikoralll9 96': Volo nteri et alj|2005l) . Recent sim- 
ulations by ISerti & Volonteri (2008 ) demonstrate that the exact 
spin distribution expected depends strongly on the nature of these 
events, and in particular how the angular momentum of the accreted 
material/coalesced black holes relate to one another. If the accretion 
is primarily 'ordered', i.e. the accreted material has angular mo- 
mentum in the same sense as that of the black hole, the black holes 
are spun up and will preferentially display high spins. If instead the 
accretion is chaotic, i.e. the angular momenta of the accreted ma- 
terial and the AGN are not typically aligned, the opposite will be 
true, leading to a strong preference for non-rotating black holes. 

Measuring black hole spin requires knowledge of the radius of 
the innermost stable circular orbit (ISCO). Under the assumption 



that the inner accretion disc extends in to this radius, the ISCO can 
be measured through study of the reflected emission produced as 
the optically thick disc is irradiated by the hard X-rays generated 
in the corona, most likely via Compton up-scattering. This emis- 
sion contains a combination of backscattered continuum emission 
and atomic features lines, the most prominent of which is often 
the iron Kq emission line (6.4—6.97 keV, depending on ionisation 
stat e) owing to its high cosm ic abundance and fluorescent yield (see 
e.g. lGeorge&Fabianlll99lh . The reflected emission from the in- 
ner disc is blurred and broadened by the strong Doppler and rel- 
ativis tic effects inherent to material in close orbit around a black 
hole iFabian er al. 1989; Laor 19913. This can results in broad and 
skewed observed emission lines profiles, even though the atomic 
emission is intrinsically narrow, and if the amount of relativistic 
blurring can be constrained, the inner radius of the disc can be 
quantified. Active galaxies do often display broad emission fea- 
tures that ca n be associated with relativistically broadened iron 
lines (see e.g. Tanaka et alj 19951 Nandra et alj|l99sj . iFabian et al.l 



I2OOOI iNandra et alj|2007l IFabian et alJ^OOglT and although ques- 
tions have been raised over whether this is the correct interpreta- 
tion for these features, wit h absorption ba sed models propose d as 
an alternative jMiller et all2008. . .2009l) . in lWalton et al] 1 I2OI2I) we 
present a strong argument in favour of the relativistic disc reflec- 
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tion origin througli comparison witli tlie similar features frequently 
observed in Galactic black hole binaries. 

Although the iron Ka line is frequently the most prominent 
reflected feature, for a broad range of ionisation states the full re- 
processed emission spectrum also contains a host of other emission 
lines from lighter elements at lower energies. In addition, reflection 
spectra also display a broad emission feature peaking at ^^30 keV, 
referred to as the Compton hump, due to the combined effects pho- 
toelectric absorption of photons at lower energies and Compton 
down-scattering of photons at higher energies. In order to self con- 
sistently treat the reflected emission arising from the inner regions 
of the accretion disc, the same processes that broaden and skew the 
iron line must be applied to all of these features. These effects can 
blend the soft complex of reflected emission lines together into a 
smooth emission feature, offering a natural explanat ion for the soft 
excess frequently observed in the spectra of AGN jCrummv et al.l 
I2OO6I) . Therefore, where possible, the broadband spectra covering 
all the key reflected emission features should be considered when 
attempting to constrain the spin of AGN, and in principle it should 
be possible to constrain the spin from the broadband spectrum even 
for sources that lack prominent broad iron lines, assuming the soft 
excess is indeed associated with reflection. 

The best sources with which to perform such analysis and rig- 
orously test the reflection interpretation are therefore 'bare' active 
galaxies, sources with little or no intrinsic absorption to compli- 
cate the determination of the r eflected emission. No t able examples 
of such sources are Fairall9 ( Gondoin et alj 1200 ll : ISchmoll et al. 



2004 



20091 : lEmmanoulopotilos et all l20Ilh . A rk 120 iVaughan et al 



Nardini etal J l201lh . Ton SI 80 jVaughan et al.l I2OO2I : 



Nardini et al .I I20I2I). and the extr eme narrow line Seyfert 1 galaxy 



IH 0707-495 l lFabian et alj|2009l) . 

However, even for apparently bare AGN, interpretations other 
than reflection have been proposed to explain the soft excess. Al- 
though bare AGN show little or no absorption by partially ionised 
material when studied by high resolution grating spectrometers, it 
was suggested that the broad, overall curvature observed at low en- 
ergies could still be caused by such absorption in an extremely tur- 
bulent outflow, in which large velocity broadening smears out the 
discrete atomic features jGierlihski & Don3l2004l : wliddleton et al.l 
l2007t) . However, more detailed simulations show that outflows 
driven from the disc are not actually able to p r ovide the ve- 
locity broadening required jSchurch & DondboOTl : [Schurch et al.l 
120091) . More recently, complex hybrid Comptonisation interpreta- 
tions have successfully been applied to the data, in which the soft 
excess arises due to Comptonisation by thermal electrons within 
the inner disc, possibly existing as a skin-like layer on the surface 
of the disc, while the hi gher energy continuum arises in a separate, 
more traditional corona l lDone et al.ll2012l : |jin et al.l20I2h . 

Nevertheless, there is very good evidence to suggest that 
the soft excess does indeed arise due to reprocessing of the pri- 
mary continuum. First of all, the energy of the soft excess is ob- 
served to be remarkably constant over many orders of magnitude 
in black hole m ass and al so over a fair range in Eddington ratio 
jGierlinski & Do ne 2004!: iMiniutti et al.l [2009b , which would re- 
quire a roughly constant electron temperature in the Comptonisa- 
tion interpretatio n. Indeed, roughly con stant electron te mperatures 
were obtained bv lDone et al. I ( I2OI2I) and ljin et alj ( l2012h . The tem- 
perature of the accretion disc is dependent on both black hole mass 
and Eddington ratio, and should vary substantially over the range 
of black hole masses and luminosities observed. It is probably not 
unreasonable to expect that some variation should also be seen in 
the temperature of the disc electrons over these ranges. The con- 



stant energy of this feature therefore strongly hints at an atomic 
origin which, given the inability of absorption to reproduce the soft 
excess highlighted previously, supports a reflection interpretation. 

In addition, one of the key observational discoveries over 
the last couple of years is the detection of lags between the 
soft and hard X-ray emission, which demonstrate that, on short 
timescales, the sof t excess responds t o changes in the high 
energy continuum 1 Zoghbi et al.l 20 id: Zoghbi & FabiarJ 1 20 III : 



lEmmanouloDoulos et al] 2011 ; de Marco et alj|201 ll) . as predicted 
by the reflection based interpretation. Although there have been 
claims that absorption based interpretations can also reproduce 
these lags via reflection fro m the more distan t obscuring clouds 
in the case of IH 0707-495 jMiller et alj|2010h . this requires both 
a special source geometry and viewing alignment, while evidence 
for such lags have no w been seen in a large number of sources 
jPe Marco et alj20I2l) . arguing strongly against this interpretation. 
Crucially, a similar lag has also now been detected between the high 
energy continu um and the broad c omponent of the iron line profile 
in NGC4I5I jZoghbi et al.ll2012h . strongly suggesting that such 
lags are indeed characteristic for reflected emission. These lags also 
appear to argue strongly against the Comptonisation interpretation, 
as fluctuations in the corona should not precede fluctuations in the 
inner disc. 

In this paper, we present a reflection based analysis of a sam- 
ple of bare AGN, with the purpose of both testing the robustness 
of the reflection interpretation, and where possible providing initial 
spin constraints for the AGN considered. We specifically focus on 
Suzaku observations in this work in order to utilise its combination 
of broadband spectral coverage and large collecting area. The pa- 
per is structured as follows: section |2] describes our data reduction 
and sample selection procedure, and section [3] details our spectral 
analysis. We discuss our results in section|4]and finally summaries 
our conclusions in section |5] 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Data Reduction 

We reduced all the Suzaku observations of type 1 AGN publicly 
available prior to October 2010. Using the latest HEASOFT software 
package we processed the unfiltered event files for each of the XIS 
CCDs and editing modes operational in the respective observations, 
following the Suzaku Data Reduction Guid^B We started by creat- 
ing new cleaned event files for all the XIS detectors operational 
during each observation (note that XIS2, one of the front illumi- 
nated detectors, experienced a charge leak in November 2006 and 
has not been in operation since) by re-running the Suzaku pipeline 
with the latest calibration, as well as the associated screening crite- 
ria files. In addition to the standard XIS event selection criteria, we 
also required that the cut-off rigidity for clean events be greater 
than 6GeV. The exact size of the circular source regions used 
varied from source to source, depending on e.g. the presence of 
other nearby sources, but were typically ~3.5' in radius, and back- 
ground regions were selected from the surrounding areas on the 
CCD free of any contaminating sources, with care taken to avoid 
the calibration sources in the comers. XSELECT was used to extract 
spectral products from these event files, and responses were gener- 
ated for each individual spectrum using the XISRESP script with a 
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Table 1. Observational details for the sample of 'bare' AGN compiled in this work, ordered in terms of the total XIS counts recorded. The uncertainties quoted 
on the XIS countrates ar e the 1 a uncertainties , while the PIN source rates are quoted as a percentage of the total rate recorded by the detector. Galactic column 
densities are taken from lKalberla et al] fcOOSi) , and positions and redshifts are compiled from the NASA Extragalactic Database (NED). 



Source 


RA 


DEC 




-'^H,Gal 


OBSID 


UUa. UdLc 


Exposure 


FT VTS Rntp 
r^i .^1.3 Rate 


PIN % 




(h:m:s) 


(d:m:s) 




(lO^O cm-2) 






(XIS/PIN; ks) 


(ct/s) 




MrkSuy 


20:44:09.7 


-10:43:25 


0.0344 


4.25 


701093010 


25/04/2006 


22/15 


3.44b ± 0.008 


30 












701093020 


14/10/2006 


24/22 


4.416 ± 0.008 


25 












701093030 


15/11/2006 


18/17 


4.42 ±0.01 


24 












701093040 


27/11/2006 


28/28 


3.821 ±0.008 


23 


3C 382 


i O.J J .yjj .t 


-t-32'41 '47 


U.UJ / o / 


6 98 




01/11 /9007 


120/1 14 


Q 07^^ -t- 004 


22 


Mrk 335 


00:06:19.5 


-1-20:12:10 


0.0258 


3.56 


701031010 


21/06/2006 


151/132 


1.569 ±0.002 


7 


Fairall 9 


01:23:45.8 


-58:48:20 


0.0470 


3.16 


702043010 


07/06/2007 


145/127 


1.963 ±0.003 


15 


IH 0419-577 


04:26:00.8 


-57:12:00 


0.1040 


1.26 


702041010 


25/07/2007 


179/161 


1.453 ±0.002 


14 












ir\Ar\^ Ar\-i A 

7040640 10 


1 0/0 1/20 10 


1 1 r\ /I ni 

1 19/10/ 


0.953 ± 0.002 


7 


Ark 564 


22:42:39.3 


-1-29:43:31 


0.0247 


5.34 


7021 17010 


26/06/2007 


89/1 10 


2.678 ± 0.004 


7 


Ark 120 


05:16:11.4 


-00:08:59 


0.0327 


9.78 


702014010 


01/04/2007 


91/104 


2.240 ± 0.004 


16 


3C 390.3 


18:42:09.0 


-1-79:46:17 


0.0561 


3.47 


701060010 


14/12/2006 


91/92 


2.083 ± 0.004 


22 


PKS 0558-504 


05:59:47.4 


-50:26:52 


0.1372 


3.36 


701011010 


17/01/2007 


21/18 


1.415 ±0.006 


8 












-7m All r\"ir\ 

70101 1020 


18/01/200/ 


19/17 


2.001 ± 0.007 


6 












701011030 


19/01/2007 


21/18 


1.251 ± 0.006 


4 












701011040 


20/01/2007 


20/16 


2.049 ±0.007 


7 












701011050 


21/01/2007 


20/15 


2.078 ± 0.007 


5 


NGC 7469 


23:03:15.6 


+08:52:26 


0.0163 


4.45 


703028010 


24/06/2008 


92/85 


1.258 ±0.003 


19 


Mrk 110 


09:25:12.9 


-f-52:17:ll 


0.0353 


1.30 


702124010 


02/11/2007 


91/80 


1.414 ±0.003 


13 


SwiftJ0501.9-3239 


05:19:35.8 


-32:39:28 


0.0124 


1.75 


703014010 


1 1/04/2008 


36/34 


2.362 ± 0.006 


22 


Mrk 841 


15:04:01.2 


+10:26:16 


0.0364 


2.22 


701084010 


22/01/2007 


45/41 


0.950 ± 0.003 


12 












701084020 


23/07/2007 


45/60 


0.961 ±0.003 


12 


Ton SI 80 


00:57:19.9 


-22:22:59 


0.0620 


1.36 


701021010 


09/12/2006 


121/102 


0.820 ± 0.002 


4 


PDS 456 


17:28:19.8 


-14:15:56 


0.1840 


19.6 


701056010 


24/02/2007 


174/209 


0.381 ±0.001 


2 


1H0323-H342 


03:24:41.1 


+34:10:46 


0.0610 


12.7 


704034010 


immm^ 


76/69 


0.559 ±0.002 


11 


UGC6728 


11:45:16.0 


+79:40:53 


0.0065 


4.42 


704029010 


06/06/2009 


46/39 


0.897 ±0.003 


14 


Mrk 359 


01:27:32.5 


+19:10:44 


0.0174 


4.26 


701082010 


06/02/2007 


98/96 


0.307 ±0.001 


4 


MCG-2-14-9 


05:16:21.2 


-10:33:41 


0.0285 


7.13 


703060010 


28/08/2008 


131/120 


0.217 ±0.001 


5 


ESO548-G081 


03:42:03.7 


-21:14:40 


0.0145 


2.29 


704026010 


03/08/2009 


39/33 


0.661 ±0.003 


15 


Mrk 1018 


02:06:16.0 


-00:17:29 


0.0424 


2.43 


704044010 


03/07/2009 


37/35 


0.548 ± 0.003 


13 


RBS 1124 


12:31:36.4 


+70:44:14 


0.2080 


1.52 


702114010 


14/04/2007 


79/83 


0.241 ±0.001 


10 


IRAS 13224-3809 


13:12:19.4 


-38:24:53 


0.0658 


5.34 


701003010 


26/01/2007 


166/159 


0.102 ±0.001 




IH 0707-495 


07:08:41.5 


-49:33:06 


0.0406 


4.31 


700008010 


03/12/2005 


121/137 


0.060 ±0.001 




IRAS 05262-1-4432 


05:29:55.5 


+44:34:39 


0.0322 


32.1 


703019010 


12/09/2008 


68/66 


0.088 ±0.001 





medium resolution. The spectra and response files for all the front- 
illuminated instruments, XISO, XIS2 (where operational) and XIS3 
were combined using the FTOOL ADDASCASPEC. Since we are in- 
terested in the average spectral properties of active galaxies in this 
work, where a source was observed on more than one occasion the 
individual spectra obtained from these observations were also com- 
bined into a single, averaged spectrum; for these sources, the spec- 
tra obtained from each of the individual observations were found to 
be broadly similar in each case. Finally, we grouped the spectra to 
have a minimum signal-to-noise (S/N) of 5 per energy bin with the 
SPECGROUP task (part of the XMM-Newton SAS), to allow the use 
of minimization during spectral fitting. During modelling, we 
do not consider the XIS spectra below 0.6 keV and also exclude the 
1.7-2.5 keV energy range owing to calibration uncertainties, and 
we primarily make use of the front illuminated spectra, owing to 



its superior performance at high energies, although the back illumi- 
nated spectra are routinely checked for consistency. 

For the HXD PIN detector we again reprocessed the unfiltered 
event files for each of the observations considered following the 
data reduction guide. Since the HXD is a coUimating instrument 
rather than an imaging spectrometer, estimating the background re- 
quires individual consideration of the non X-ray instrumental back- 
ground (NXB) and cosmic X-ray background (CXB). The response 
and NXB files were downloaded for each observatior0, for all of 
the observations considered here the higher quality 'tuned' (Model 
D) background was available. Spectral products were generated us- 
ing the HXDPINXBPI script provided by the Suzaku team, which 



^ http://www.astro.isas.ac.jp/suzaku/analysis/hxd/ 
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includes a si mulated contribution from the CXB using the form of 
iBoldll ( Il987l) . As with the XIS detectors, we combined the indi- 
vidual spectra of sources with multiple observations into a single, 
averaged spectrum, and rebinned the data so that each energy bin 
had a minimum S/N of 3, lower than the minimum S/N required 
in the XIS spectra owing to the much lower PIN source rates, but 
still sufficient to allow the use of minimization. Nominally we 
consider the PIN spectrum over the full 15-70 keV energy band, 
but in reality the exact energy range considered varies from source 
to source, depending on the highest energy at which each source is 
reliably detected. 



2.2 Sample Selection 

In this work, we are interested in 'bare' AGN, i.e. active galax- 
ies with little or no intrinsic absorption, that display soft excesses. 
In order to identify such sources, we modelled the --^2.5-10 keV 
energy range with a powerlaw continuum modified by Galactic 
absorption (see Table [T] f or the Galactic neutr al column densities 
adopted, as determined bv'K alberla et alJlOOSi) . modelled with the 
TBABS absorption code l lWilms et alj|200(ll : the appropriate solar 
abundances for this model are used throughout); the iron K band 
(generally taken here as 4-7 keV) was excluded to prevent any iron 
emission present in the spectrum from influencing our continuum 
determination. We then extrapolated this continuum back through 
the 0.6-2.0 keV data and selected sources that displayed a clear, 
and relatively smooth {i.e. not domina ted by diffuse ther mal emis- 
sion, as is the case for e.g. NGC 1365: IWang et alj2009l) excess of 
emission over that predicted at lower energies. The resulting sample 
contains 25 sources that met our selection criteria, the observational 
details of these sources analysed are listed in Table[T] Note that this 
is largely based on a visual inspection of the spectra, so some of the 
sources selected may (and indeed do) display some mild, partially 
ionised absorption, but even in this minority of cases this absorp- 
tion does not have a large effect on the observed spectrum. 



3 SPECTRAL ANALYSIS 

Here we detail our spectral analysis of the sample outlined in the 
previous section (see Table [T). In section 13.11 we outline the ba- 
sic approach we adopt for all the sources, and in section IT2l we 
provide additional, detailed information on a source-by-so urce ba- 
sis. S pectral modelling is performed with XSPEC vl2.6.0f jAmaudl 
Il996h . and throughout this work the uncertainties quoted on model 
parameters are the 90 per cent confidence limits for one parameter 
of interest, unless stated otherwise. We also adopt as standard the 
cross-normalisation constant of either Cpin /xiso = 1.16 or 1.18 be- 
tween the XISO and PIN spectra (the former for observations with 
the nominal pointing position centred on the XIS detectors, and the 
latter on the HXD detector), as recommended by the HXD calibra- 
tion tearrffl However, since we are using combined FI XIS spectra 
in this work (XISO, XIS2 and XIS3 where available) we take the 
extra step of comparing this with the XISO spectrum on a source by 
source basis in order to calculate the correct normalisation constant 
for use between the combined XIS and PIN spectra, Cpin /xis> us- 
ing equation [T] 



http://www.astro.isas.jaxa.jp/suzaku/doc/suzakumemo/suzakumemo- 
2008-06.pdf 



CpiN/XIS — CxiSO/XIS ^ CpiN/XISO (1) 

We calculate Cxiso/xis by simultaneously modelling the com- 
bined FI XIS and XISO spectra with a phenomenological broken 
powerlaw model for the broad 0.6-10.0 keV continuum (excluding 
again the energy ranges 1.7-2.5 keV due to calibration uncertain- 
ties and 4.0-7.0 keV to exclude any iron emission), with all the 
physical parameters linked for the two spectra, and only the nor- 
malisation constant allowed to vary. The values of Cpin/xis we 
obtain are presented in Table [2] However, while this gives the ap- 
propriate normalisation constant for the situation in which the PIN 
background is known perfectly, it does not account for any sys- 
tematic uncertainties in the modelled background currently avail- 
able, which are estimated by the HXD team to be 3 per cenjj there- 
fore we also investigate the consequences of allowing the XIS/PIN 
cross-normalisation parameter to vary within the range equivalent 
to ±3 per cent of the background, these ranges and the values ob- 
tained are quoted in Table [3] 



3.1 Basic Approach 

In this work we are investigating the reflection based interpreta- 
tion for the broad band X-ray spectra of active galaxies, i.e. that 
there are two main continuum components: a powerlaw-like com- 
ponent (PLC) most likely arising via electron scattering in some 
kind of corona, and a reflected component (RDC) arising due to ir- 
radiation of the disc by the PLC, which is modified by the strong 
gravitational effects present close to the black hole. This interpre- 
tation potentially provides a natural explanation in the form of the 
RDC for the soft excesses, broad iron lines and the hard excesses 
often displayed by AGN, and forms the basis of our spectral mod- 
elling. We try to take a broadly systematic approach, as detailed 
below, while also considering the specific details of each source in- 
dividually, which are presented in section [T2l although as the data 
quality declines so in general does the level of detail to which each 
source can be considered. The spectral complexities displayed by 
each source are highlighted in Figure [T] in which we show the ra- 
tio of the XIS (black) and PIN (red) data to the absorbed powerlaw 
model discussed previously, using the cross-calibration normalisa- 
tions calculated with equation[T] 

Naturally the PLC is modelled simply with a powerlaw emis- 
sion component, with variable flux and photon index (F). To model 
the RDC component we ma ke use of the physicall y self-consistent 
reflection code REFLIONX ( IRoss & Fabiarill2005l) . This model is 
specifically designed to account for the reprocessed emission from 
the cool, partially ionised accretion discs in AGN, intrinsically in- 
cluding both the reflected continuum and the imprinted atomic fea- 
tures, in particular the iron K-shell transitions. The key param- 
eters of REFLIONX are the photon index of the ionising contin- 
uum (assumed to be powerlaw-like), and the iron abundance (Ape) 
and ionisation state of the reflecting medium. The latter is quanti- 
fied as the classic ionisation parameter for a photoionised plasma, 
^ = L/nR^ , where L is the ionising luminosity, n is the num- 
ber density of the reflecting medium and R is its distance from the 
primary X-ray source. Throughout this work, we require that the 
photon index of the ionising continuum be that of the PLC, while 



* http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/ 
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Figure 1. Data/model ratio plots of Suzaku XIS (front illuminated; black) and PIN (red) spectra to a powerlaw continuum model, as described in the text 
( section im . for the compiled sample. Each source displays a soft excess, and a number also show evidence for broad iron emission components, and/or for 
'hard' excesses above ~10 keV, although these features are by no means ubiquitous. The data shown have been re-binned for plotting purposes. 
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Figure[T] Continued; note that the photon indices of IRAS 13224-3809 and IH 0707-495 have both been fixed at 2.7 (see section [X2l . Furthermore, there are 
no robust PIN detections for either of these sources, or for IRAS 05262-H4432. 



the iron abundance and ionisation parameter are free to vary (unless 
stated otherwise, see section |T2ll. 



In order to account for the strong gravitational and Doppler ef- 
fects that the emission from material orbiting close to the black hole 
will experience, we modify this reflected emission with RELCONV 
jPauser et al.ll201(il) . one of the most sophisticated relativistic con- 
volution kernels. Under the assumption that the inner radius of the 
accretion disc is coincident with the ISCO, this kernel allows the 
black hole spin to be constrained directly from the RDC, and even 
allows for the intriguing possibility of black holes spinning in a 
retrograde sense with respect to the material orbiting in their accre- 
tion discs. The key free parameters of RELCONV are the spin of the 
black hole, the inclination of the disc (i), and the radial emissivity 
profile of the disc. We assume a single powerlaw emissivity profile 
for simplicity, i.e. e(r) oc r~'', where the emissivity index q is a 
free parameter. For a simple Newtonian 'lamp post' like geometry, 
in which the illuminating source is compact and centrally located, 
the predicted emissivity index is g = 3. However, this does not take 
into account any of the effects of general relativistic light bending, 
which focuses a larger fraction of the emission onto the inner disc 
if the X-ray source is als o located in a region of extreme gravity 
jMiniutti & Fabianll2004l) . or of relativistic time dilation. Both of 
thes e effects g enerally serve to steepen the inner emissivity pro- 
file JWilkins & Fab ian 20l"lb, therefore throughout this work we 
require that q ^ 3. The inner and outer radii of the disc are fixed 
respectively at the (spin dependent) ISCO, and at 400 -Rg (^g = 
GM/c\ the maximum value accepted by the model. It is also pos- 



sible to include the effects of limb brightenin g/darkening , and we 
assume the latter following the calculations of lLaod ( Il99 1[) . 

With this RDC we attempt to simultaneously account for the 
broad curvature across the 0.5-10.0 keV energy range that causes 
the appearance of a soft excess, any relativistically broadened 
iron emission detected and also any hard excess displayed by the 
AGN. However, this is a complex, multi-parameter model compo- 
nent, and degeneracies can arise between various parameter com- 
bin ations, particularly i n the relativistic blurring parameters (see 
e.g. iNardini et al.ll201ll) : this can be especially problematic with 
moderate quality data, or when the features associated with reflec- 
tion are weak. One of the main aims of this work is to provide 
where possible an initial, albeit model dependent constraint on the 
black hole spin. Therefore, in cases where it is not possible to in- 
dependently constrain all the blurring parameters, we fix e.g. the 
inclination to 45°and/or the emissivity index to g = 3 in an attempt 
to provide at least a loose spin constraint. This was only necessary 
for the minority of cases (7 AGN in total) even including our anal- 
ysis in which Cpin /xis was allowed to vary, although in 4 of these 
sources we were still unable to provide any constraint on the spin. 
Furthermore, in a few individual sources certain parameters might 
also be fixed for other reasons, e.g. the inclination being known 
a priori from radio constraints. Full details are provided in section 
El 

In addition to their broadband continuum emission, active 
galaxies also very frequently display evidence for narrow emis- 
sion features consistent with neutral iron (6.4 keV). Our general 
assumption is that these features also arise through reflection of the 



© 0000 RAS, MNRAS 000.[Tll2T] 



Suzaku Observations of 'Bare ' AGN 7 



primary X-ray emission, but by cold, distant material rather than 
by the accretion disc, possibly associated with the torus-like struc- 
ture that obscures the central X-ray sources of AGN viewed at high 
inclination angles. Therefore, where we see evidence for such a 
feature, our default procedure is to include a second, cold reflec- 
tion component. This was also modelled with REFLIONX, with ^ 
= 1.0 (the lowest ionisation state accepted by the model), while 
we required the iron abundance to be the same as that of the disc 
reflection component, and we again required that the photon in- 
dex of the ionising continuum to be the same as that of the PLC. 
This second reflector accounts for the narrow emission from neutral 
iron, and where present may also contribute some emission to the 
hard excess. Less common than narrow emission from neutral iron, 
which is almost observed ubiquitously, AGN also sometimes dis- 
play narrow emission features consistent with highly ionised iron 
(Fe XXV and XXVI; 6.67 and 6.97 keV respectively). Where we 
find evidence for these additional narrow features, we simply ac- 
count for them with narrow (a — 10 eV) Gaussian emission lines. 

Finally, all emission components are modified by neutral ab- 
sorption due to the Galactic ISM. Again, we make use of the TBABS 
absorption code, and fix the column density at the Galactic value 
quoted in Table [T] Given that we have deliberately selected 'bare' 
active galaxies, we do not include any contribution from neutral ab- 
sorption intrinsic to the source. However, despite our source selec- 
tion method, high resolution grating spectra of some of the sources 
included in this work do show evidence for some weak, partially 
ionised absorption intrinsic to the source. Where this is the case, 
we i nclude an ionised absorber using the XSTAR photoionisation 
code jKallman & Bautistal200ll) to account for the main features in 
our lower resolution CCD spectra. The free parameters of the ab- 
sorption model used here are the ionisation state (again quantified 
as ^), the column density and the line-of-sight velocity of the ab- 
sorbing medium; the abundances of the absorbing medium are as- 
sumed to be solar, and its turbulent velocity is fixed at 200 km ^ . 
In addition, there has been growing evidence for AGN displaying 
hig hly blueshifted absor ption features from highly ionised iron (see 
e.ie. lTombesi et al .1201(1 typically interpreted as evidence for rapid 
outflows (although see lGallo & Fabianll201ll for an alternative ex- 
planation). Where we find evidence for similar features, we also 
treat these with XSTAR, although in these cases we will use grids 
specifically designed to treat absorption by highly ionised material, 
as detailed below. 

3.2 Individual Source Details 

Here, we provide further information specific to our individual 
modelling of each source, including details on any additional emis- 
sion or absorption features included beyond the basic PLC-fRDC 
interpretation. The key reflection parameters obtained from our 
spectral modelling are quoted in Table |2] (Cpin /xis fixed at the 
values obtained from equation[Tll and in Table[3](CpiN /xis frse to 
vary within the range defined by the 3 per cent systematic uncer- 
tainty on the modelled PIN background). We also show 1-D confi- 
dence contours of the spin parameter for the full range considered 
for each source in Figure |2] based where possible on the results 
presented in Table[3] Note that these contours have been generated 
with the inclination limited to i ^ 75° for sources in which this 
parameter is free to vary. 

Mrk 509: In addition to the basic PLC-l-RDC interpretation for 
the continuum we also include a moderately ionised absorber at the 
redshift of the host galaxy, with log^ = 2.03 ± 0.04 erg cm s~^ 
and Nn= l-Oloj x 10^^ atom cm . The ionisation obtained is 



consistent with that of the most prominent absorbing zone found 
in the detailed RGS analysis presented by lOetmers et al.l 1 I2OIII) 
as part of the substantial multi-wavelengt h observing campaig n re- 
cently carried out on this source (see also lKaastra et aljbol ih . but 
the column obtained here is larger This may be due to our sim- 
plistic, single zone treatment of the absorber, but the inclusion of 
a second ionised absorbing zone does not offer a substantial im- 
provement. We also find a narrow emission feature consistent with 
neutral iron, so include a distant reflector, and additional emission 
consistent with a Fe XXV. The latter is modelled with a narrow 
Gaussian emission component, and we find an equivalent width of 
£Wxxv = 17^ 10 eV; the inclusion of this feature provides an im- 
provement of Ax^ = 23 for an additional degree of freedom. Here, 
we find that the normalisation between the XIS and PIN spectra 
plays an important role in determining the interplay between the 
PLC and RDC continuum components, and hence in determining 
the reflection parameters. When Cpin/xis is allowed to vary the 
inclination obtained is significantly larger, although still reason- 
able for an unobscured AGN, and the spin is lower and less con- 
strained. Furthermore, the ionisation of the absorber increases to 
log^ = 2.38lH^ erg cm s-\ 

3C382: As one of the few broad line radio galaxies (BLRGs) 
included in this sample, it is possible to obtain an independent con- 
straint on the inclination at which the source is viewed at from 
radio observations of its jet, assuming that the jet axis is aligned 
perpendicular to the plane of the inner disc. Although it is difficult 
to rigorously test this assumption, we note that the Galactic binary 
system XTE J1550-564, one of the few cases in which it has been 
possible to independently constrain the or bit and jet inclinations, 
appea rs to be consistent with this picture jSteiner & McClintockl 
12012b . Therefore, when applying the PLC-l-RDC continuum, we 
fixed the incl ination aii = 40° , cons istent with the radio constraint 
presented bv lGiovannini et al] ( l200lh . 

In addition to this underlying continuum, there is also a nar- 
ro w feature associated w ith neutral iron, and as initially reported 
bv lSambruna et al] ( |201 1[) . evidence for a further narrow feature at 
'~-^7.5 keV. We include a distant reflector to account for the former, 
and the latter is modelled with a narrow Gaussian emission line. 
The inclusion of this Gaussian results in an improvement of Ax^ 
= 19 for 2 additional degrees of freedom, and we find a (restframe) 
line energy of Eq = 7.53 ± 0.04 keV, and an equivalent width of 
EWg = 19 ± 9eV; the energy of this feature precludes an asso- 
ciation with iron, but is consistent with moderately ionised nickel. 
Finally, we also include a moderately ionised absorber at the red- 
shift of the galaxy, similar to that found by iTorresi et alj ^20 id) 
in their analysis of the high resolution XMM-Newton RGS spec- 
trum of this source, with log^ = 2.5 ± 0.1 erg cm s~^ and A'h 
= lAtol X 10^^ atom cm"^. As with Mrk 509, the normahsa- 
tion between the XIS and PIN spectra plays an important role in 
determining the interplay between the PLC and RDC continuum 
components in this case, and hence in constraining the reflection 
parameters. We find that when Cpin/xis is allowed to vary the 
iron abundance obtained is much lower, and the spin and the emis- 
sivity index are essentially unconstrained. Therefore in Table [3] we 
present the results obtained with the index fixed at the Newtonian 
prediction (q = 3), although the spin is still only very weakly con- 
strained. 

Mrk 335: In addition to the basic PLC-l-RDC model adopted 
in this work, we also see evidence for narrow emis sion features 
consiste nt with neutral iron and Fe XXV, similar to [Patrick et al.l 
( l201ld) . The former is modelled with a distant reflector, and the 
latter is simply modelled with narrow (a = 10 eV) Gaussian line 
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Table 2. Key Parameters obtained for the reflection based models constructed for the compiled sample (see section [X2l for details). Parameters in parentheses 
have not been allowed to vary, and where we were unable to constrain the black hole spin this is indicated with a 'U' . 



Source 


En. Range 
(keV) 


Cpj^ /xis 


r 


-4fo 
(solar) 


5 

(erg cm s^^) 


q 


i 

(°) 


a* 


xl 

(X^/d.o.f) 


Mrk:509 


0.6^4.0 


1.17 


2.04 ± 0.01 


0.5 ± 0.1 


1 -70+30 


> 7.4 


< 18 


00+0 02 
0-86Io:oi 


1314/1333 


3C382 


0.6-53.0 


1.14 


1.81 ± 0.01 


> 5.1 


5001^^0 


> 6.1 


(40) 


„ r.r--f 07 


1366/1317 


Mrk 335 


0.6-36.0 


1.17 


2.16 ± 0.01 


1.4 ± 0.2 


220 ± 10 


> 4.9 


501° 


0-83Io;i3 


1233/1152 


Fairall 9 


0.6-39.0 


1.16 


1.99 ± 0.01 


1.1 ± 0.2 




> 3.5 


.1^ + 13 
45I9 


> 0.64 


1276/1253 


IH 0419-577 


0.6^8.0 


1.19 




0.9 ± 0.1 




_ .+0 2 


r--, 4-4 

Slle 


> 0.88 


1384/1294 


Ark 564 


0.6-22.0 


1.20 


2.52 ± 0.01 


1 


540l*° 


> 6.2 


^ .-1-1 


^ ^-^4-0 01 


1081/1011 


Ark 120 


0.6-53.0 


1.22 


2.13 ± 0.01 


2.7 ± 0.6 


9Ii 


„ 0+2 7 
7-012:2 




o-8iio:i8 


1316/1177 


3C 390.3 




1 1 ^ 
i.io 


l.OD ± U.Ul 


n 1+1.2 
3-1-0.5 


840_220 


(3) 


(35) 


u 


1 /I T ^ n 


PKS 0558-504 


0.6-24.0 


1.13 




r. r.4-0 3 

0-9Io.i 




4.0 ± 0.7 


(45) 


> 0.80 


990/1022 


NGC 7469 


0.6-52.0 


1.19 


^ .4-0 03 
1 S4 

-^•°^-0.02 




2001^0 


> 4.6 


< 54 


> 0.96 


1262/1139 


Mrk 110 


0.6^5.0 


1.19 


1 nc+0.02 


0.7 ± 0.2 


oi n+220 
310_80 


> 7.4 


-91 +4 
31-6 


> 0.99 


1184/1115 


Swift J050 1.9-3239 


0.6-36.0 


1.22 


^•^'^-0.03 


i-8to;5 


200l^^ 


> 5.1 


< 48 


> 0.96 


1025/1056 


Mrk 841 


0.6-53.0 


1.19 


-X.^^ ^'^ 

i-«o_o.oi 


1.0 ±0.2 


-4-90 

210I70 


A 1 +^-" 


45I5 


> 0.56 


1089/1053 


Ton SI 80 


0.6-23.0 


1.16 


2.36 ± 0.01 


n n + 0.2 

0-9-0.1 


2801^^1 


> 8.1 


60^1 


„ +0,02 
O-9II0.09 


876/838 


PDS 456 


0.6-17.0 


1.15 


r, or,+0.03 

2-30Io,oi 


> 8.4 


59I4 


c 0+1.8 

5-9-1.5 


70lg 


> 0.97 


829/826 


IH 0323-1-342 


0.6^2.0 


1.25 


-'-■^-'--0.01 


0.8 ± 0.2 


250l« 


(3) 


(45) 


> 0.48 


864/922 


UGC 6728 


0.6-26.0 


1.27 




'-0.3 


190l??o 


6-81?:! 


< 55 


> 0.95 


877/885 


Mrk 359 


0.6-21.0 


1.15 


i-89io:M 


l-5lo.6 


+32 

21-16 


> 4.1 


47 ±6 


o-66lE!:^g 


820/833 


MCG-2-14-9 


0.6-37.0 


1.19 


1.89 ±0.02 


(1) 


< 10 


(3) 


(45) 


U 


802/804 


ESO548-G081 


0.6-36.0 


1.23 


1.70 ±0.03 


3 5+4-1 


C.70+56O 
3'U_3go 


(3) 


(45) 


U 


853/845 


Mrk 1018 


0.6^1.0 


1.21 


1 94+0-04 
^•='^-0.03 


2 0+^''' 
^■'J-O.T 


c + lO 

0_4 


> 3.9 


4= + 14 


57+0-31 


681/721 


RBS 1124 


0.6-23.0 


1.22 


J^-»O_0.02 


2 9+l'5 


511^ 


> 8.4 




> 0.98 


661/668 


IRAS 13224-3809 


0.6-7.6 




(2.7) 


(20) 


22 ± 3 


°-l-0.6 


(64) 


> 0.995 


447/412 


IH 0707-495 


0.6-6.7 




(2.7) 


(10) 


531^ 


'■"-0.3 


(58) 


> 0.994 


278/236 


IRAS 05262-H4432 


0.6-7.8 




9 10+0.13 
^•^"-0.06 


(1) 


< 51 


(3) 


(45) 


U 


234/231 



profile. The inclusion of the latter feature provides a strong im- 
provement of Ax^ = 33 for 1 extra degree of freedom, although it 
does not significantly effect the parameters obtained for the main 
PLC and RDC continuum components, and we find an equivalent 
width of EWxxv ~ 29 ± 9 eV. The spin obta ined here is very sim- 
ilar to that presented bv lPatricketinhoilbl) for the same dataset, 
and is consi stent with bo t h the spin inferred from the inner radius 
obtained by iGrupe et"al] 1 I2OO8I) in their analysis of the low flux 
XMM-Newton datasets, and the spin measure d from the recent i n- 
termediate flux XMM-Newton observations bv lGallo etan ( l2012h . 

Fairall 9: In addition to the standard PLC+RDC continuum 
considered here, we also detect strong, narrow emission features 
consistent Fe I and and Fe XXVI, and a weaker emission feature 
consistent with Fe XXV. A distant reflector is included to model the 
former, and the latter two are modelled with narrow (cr = 10 eV) 
Gaussian line profiles. The equivalent obtained are SWxxv = 
13l^ eV and BWxxvi = 35l? eV for the ionised emission lines, 
and their inclusion results in respective improvements of Ax^ = 
14 and 52 for one extra degree of freedom in each case. Although 
the inclusion of the Fe XXVI feature does have a significant effect 
on some of the main reflection parameters, the best fit model ex- 
cluding this line resulted in an estimate for the inclination of 87°, 



which is unphysically high for an unobscured AGN. Conversely, 
the inclusion of the Fe XXV line does not significantly modify the 
key parameters further The spin obtai ned here is larger, alth ough 
still consistent with those obtained by Sch moll et al.l 1 I2OO9I) with 
the same dataset, by£atricket al. (2011a) with a later Suzaku ob- 
ser vation not included here (perfo rmed after October 2010), and 
by lEmmanoulopoulos et alj ^201 ih with a lo ng XMM-Newton ob- 
servation, but is lower than that presented in lPatrick et al.l ( 1201 Ibh 
again with the same dataset. 

IH 0419-577: In addition to the PLC+RDC components, we 
also detect a weak narrow emission feature consistent with neu- 
tral iron, and therefore include a second, distant reflector The spin 
constraint obtained here is consistent with that inferred both i n our 
previous work on the first Suza ku dataset dWalton et alj|2010h . and 
in the XMM-Newton analysis of lFabian et alj ^200^ 

Ark 564: The basic PLC+RDC interpretation provides a good 
fit to the data. No narrow iron emission or absorption lines are de- 
tected. The spin obtained h ere is consistent with th at infeiTed from 
the inner radius obtained bv lDewangan et af] ( l2007l) . who modelled 
the long 2005 XMM-Newton observation with blurred reflection. 

Ark 120: In addition to the basic PLC+RDC continuum, we 
also detect a narrow emission feature consistent with neutral iron, 
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Table 3. Key Parameters obtained for the reflection based models constructed for the compiled sample when Cpj^ /xis '° within the indicated 

range. Again, parameters in parentheses have not been allowed to vary, and where we were unable to constrain the black hole spin this is indicated with a 'U'. 



Source C'pin/xis (range) T Apc ^ q i a* xl 









(solar) 


(erg cm s ^ ) 




(°) 




(X^/d.o.f.) 


Mrk 509 


< 1.06 (1.04-1.30) 


1.99 ± 0.01 


0.5 ± 0.1 


-1 n + 8 

13l« 


> 4.3 


rn + 5 


0-36i;j;^'; 


1302/1332 


3C382 


< 1.03 (0.99-1.29) 


1.82 ± 0.01 


l-7Io.3 


or>r»+10 

200lyQ 


(3) 


(40) 


< 0.81 


1358/1317 


Mrk 335 


1.19 ± 0.11 (0.79-1.55) 


2.16 ± 0.01 


1.4 ± 0.2 


220 di 10 


> 4.9 


co+S 

50„9 


n 0-3 + 0.09 

0-83Io.i3 


1233/1151 


Fairall 9 


1.06 ± 0.06 (0.95-1.37) 


2.00 ± 0.01 


^ ^4-0 2 

o-9io:i 


13013^ 


_ f^+4 5 




„ 09 

0.82+^;™ 


1269/1252 


IH 0419-577 


1.27j;Q;Q5 (0.93-1.45) 


1 0-7+0.03 

l-97Io.oi 


0.8 ±0.1 


1 <in+10 
120+20 


r- n + 0.3 

5-2Ii.3 


An+7 
491 4 


> 0.89 


1380/1293 


Ark 564 


l-40to (0.79-1.61) 


2.53 ± 0.01 


-. ^-t-O 2 
J^-U-O.l 


550l*[J 


> 6.0 


^ .4-1 


0.961°;°^ 


1078/1010 


Ark 120 


< 1.04 (1.00-1.43) 


2.14 ± 0.01 


1 O+0.2 

1-8-0.3 


e + 5 
6-2 


c 1 +2.8 

5-1-1.1 


47I2 


n c.1+0.19 

0-64Io.ii 


1274/1176 


3C 390.3 


1.2ll[5;o4 (1.01-1.31) 


1.65 ± 0.01 


A t; + 1.7 


1060+350 


(3) 


(35) 


U 


1301/1258 


PKS 0558-504 


> 1.23 (0.70-1.56) 


^ .^^^0 09 


^^■'^-0.2 




4.1 ± 0.7 


(45) 


> 0.95 


985/1021 


NCjC 7469 


< 1.13 (1.01-1.37) 


1 O1+0.02 
1-»1_0.01 


1 1+0.5 


2001 


> 5.5 


(45) 


n c/i+0.13 
0-64Io.20 


1261/1139 


Mrk 1 10 


< 1.07 (0.96-1.42) 


1-921^.02 


O r-r+0.2 

0-7Io.i 


^„„ + 230 
25OI40 


> 5.4 


< 43 


n r.A> + 0.03 

0-96t^.^? 


1179/1114 


Swift J050 1 .9-3239 


> 1.25 (1.07-1.37) 


2.09 ± 0.04 


o 1+1.4 
3-1-1.2 


nnn+10 
2OOI7Q 


n n+0.4 
9-0+2.7 


09+10 

•^^-28 


> 0.99 


1019/1055 


Mrk 841 


1.14+QjQ (0.90-1.48) 


1.86 ± 0.02 


n O+0.3 
0-9-0.2 


1 nn+30 
190-110 


A 1 +2.8 
4-1-1.9 


46+5 


> 0.52 


1088/1052 


Ton SI 80 


> 1.46 (0.42-1.90) 


2.36 ± 0.01 






> 8.3 


60ti 


„ „^+0 03 

0.92+o;"4 


869/837 


PDS 456 


1.34:^0,11 (0.00-2.97) 


9 on+0-03 
^•■JlJ-O.Ol 


> 8.5 


59tf 


5-9l?:i 


70^3 


> 0.96 


829/825 


1H0323-H342 


> 1.26 (0.89-1.61) 


1.92lH^ 


0.8«-i 




(3) 


(45) 


> 0.37 


861/921 


UGC 6728 


< 1.41 (1.02-1.52) 


1.98lSj;i^ 


0.8iO-j 


< 500 


> 4.5 


(45) 


> 0.71 


877/884 


Mrk 359 


0.78lo 23 (0.33-1.97) 


1.90 ±0.03 


i-2tS:I 


91+32 


> 4.1 


47^10 




814/832 


MCG-2-14-9 


l-lSi'hM (0.54-1.85) 


1.89 ±0.02 


(1) 


< 26 


(3) 


(45) 


u 


802/803 


ESO548-G081 


< 1.21 (0.98-1.48) 


-I 7-1+0.05 
^•'^-0.02 


9 9 + 1.8 
^•^-1.1 


270^520 


(3) 


(45) 


U 


850/844 


Mrk 1018 


> 1.27(0.94-1.48) 


1 qs+003 


9 7+1.7 
^-'-1.2 


7+18 
'-5 


> 3.4 


46 ±9 




676/720 


RBS 1124 


> 1.46 (0.84-1.60) 


1.86 ±0.05 


9 7+1.8 
^-'-0.9 


OU_4o 


> 6.9 


64iL 


> 0.97 


649/667 



and hence include a second, distant reflector. The spin obtained 
here is broadly si milar to, although slig htly better constrained than, 
that obtained by iNardini et al.l ( I2OI ih with the emi ssivity index 
fixed at q = 5, but is lower than that obtained by IPatrick et al.l 
( l2011bl) with the same dataset. 

3C 390.3: Similar to 3C 382, 3C 390.3 is a BLRG, and hence 
we fix the inclination at i = 35° when applying the PLC+RDC 
continuum, consistent with the radio constraint on the inclination 
presented bv lGiovannini et al] bOOll) . In addition to this underlying 
continuum, we find evidence for a strong narrow emission feature 
consistent with neutral iron, and hence include a distant reflector 
However, owing to the weak features associated with reflection, we 
are unable to reliably constrain all the parameters of interest, and 
even when fixing the emissivity index at o = 3 we find the spin 
to be unconstrained. Note that similar to Sambruna et alj ^20091) 
we also include an unresolved Gaussian component at ~1.4keV to 
account for a feature in the background spectrum. 

PKS 0558-504: Here the basic PLC+RDC interpretation pro- 
vides a good fit to the data; no statistically compelling narrow iron 
emission lines are detected. However, when free we find that the 
inclination obtained is very high (i > 75°), unphysically high for 
an unobscured active galaxy. This is most likely due to the complex 
parameter degeneracies th at can arise when us ing multi-parameter 
reflection models (see e.g. iNardini et al.ll201lh . Therefore, we in- 
stead present the results obtained with the inclination fixed at 45°, 



similar to the rough estimation of lGliozzi et alj ( 1201 Oh based on the 
radio morphology of the source. 

NGC 7469: In addition to the PLC+RDC continuum, there is 
a strong narrow emission line consistent with neutral iron, so we 
again include a distant cold reflector We also include a weak, mod- 
erately ionised absorber at the redshift of the galaxy, based on the 
most p rominent absorption component detected by iBlustin et al.l 
( I2OO7I) when modelling the high resolution XMM-Newton RGS 
spectrum. As the absorber has an extremely weak effect on the 
spectrum in this case, we only allow the column to vary, fix- 
ing the ionisation parameter at log^ = 2.73, and find TVh = 
1.7lo'g ^ lO^'^ atom cm^'^, consistent with that obtained previ- 
ously. In this case, when allowing Cpm /xis to vary it is not pos- 
sible to constrain all the reflection parameters of interest, so we 
present the results with the inclination fixed at 45°. In addition to 
the golbal minimum presented in Table [3] there is a further, local 
minimum with a very high spin (a* > 0.98) which also satisfies 
Ax^ < 2.71 (see Fig. [3}, equivalent to the best fit obtained with 
CpiN /XIS fixed presented in Table|2] 

Mrk 110: In addition to the PLC+RDC continuum, we also de- 
tect a strong narrow emission feature consistent with neutral iron, 
and a possible narrow absorption feature at ~ 7.5 keV. We in- 
clude a distant reflector to account for the former, and model the 
latter with a highly ionised photoionised plasma using the XSTAR 
code. The absorption zone here is assumed to have solar abun- 
dances and a turbulent velocity of --^1000 km s^^, while the col- 
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umn, ionisation parameter and outflow velocity are free to vary, and 
we find Nn = 9.6t?;2 x 10^^ atom cm"^ log^ = 4.2 ± 0.3 and 
VoMt = 0.079 ± 0.005c. The inclusion of this absorber provides an 
improvement of A^^ = 21 for 3 extra degrees of freedom. 

Swift J0501. 9-3239: In addition to the basic PLC+RDC con- 
tinuum, we also detect a narrow emission feature consistent with 
neutral iron, and include a second, distant reflector. 

Mrk 841: In addition to the PLC+RDC components, we also 
detect a narrow emission feature consistent with neutral iron and 
therefore include a contribution from a second, distant reflector. 
We also include a moderately ionised absorber at the redshift of 
the galaxy similar to the most prominent absorption component de- 
tected by Longinotti et al. (2010) in their analysis of the high res- 
olution XMM-Newton RGS spectrum of this source, with log ^ = 
2.2 ± 0.2 erg cm s'^ and A^h = 2.lto.8 x 10^^ atom cm"^ Al- 
though we present the prograde spin solution in Tables |2] and [3] 
there is also a local minimum with retrograde spin (a* < —0.4) 
also just about satisfies Ax^ < 2.71 (see Fig.|2l(. Th e prograde so- 
lution is consistent with the inner radii obtained by IPetrucci et al.l 
( |2007|) in their analysis of the available XMM-Newton data. 

Ton S180: The basic PLC-l-RDC continuum provides a good 
fit to the data, no narrow emission or absorption lines are detected 
over the energy range considered. The spin constraint ob tained here 
is very similar to that presented bv lNardini etal] ilOlA . who kept 
the inclination fixed at 45° . 

PDS456: In addition to the PLC+RDC continuum, there 
is also well documented evid ence for absorption due to an ex- 
tremely high velocity outflow l lReeves et ai]|2009l) . Owing to the 
high iron abundance obtained, we follow our original approach 
and model the outflow with a photoionised plasma with a vari- 
able iron abundance, fixed to that of the reflector for consistency. 
The other abundances are assumed to be solar, and the turbulent 
veloci ty is assumed to be ~3000 km s^^ following iReeves et al.l 
( l2009h . The column, ionisation parameter and outflow velocity 
are allowed to vary, and we find A^H = 7tl X 10^^ atom cm 
and log^ = 3.521q q3, consistent with our prev i ous w ork, and 
Vout = .31 ± 0.01c, consist ent with lReeves et alj j2009l) . As dis- 
cussed in IWalton etalj l l2010l) . the high iron abundance and in par- 
ticular the high inclination obtained here are driven by the associa- 
tion of the emission feature at ~0.8 keV with reflected iron-L emis- 
sion, as opposed to arising due to complex parameter degeneracies 
as is likely to be the case for PKS 0558-504, so in this instance we 
consider the values obtained to be robust. As with Mrk 841, there 
is also a local minimum with a retrograde spin (a* < —0.5) that 
satisfies Ax^ < 2.71 (see Fig.O. However, the prograde solution 
obtained is consistent with our previous work on this dataset. Note 
finally that in this case, in order to include a PIN detection it was 
necessary to relax our rebinning requirement to a S/N of 2 instead 
of 3. 

IH 0323+342: In this case the basic PLC+RDC interpretation 
provides a good fit to the data; no statistically compelling narrow 
iron emission lines are detected. However, similar to the earlier case 
of PKS 0558-504, when free we find that the inclination obtained is 
unphysically high for an unobscured active galaxy (i = 82 ± 3°). 
Therefore, we again present the results obtained with the inclina- 
tion fixed at 45° . However, we are still not able to constrain all the 
reflection parameters of interest, so we present the results obtained 
assuming g = 3. Note that when Cpin /xis is free to vary, in ad- 
dition to the prograde solution presented in Table [3] there is also 
a local minimum with a retrograde spin (a* < —0.4) which also 
satisfies Ax^ < 2.71 (see Fig.|2l(, similar to Mrk 841. 

UGC6728: In addition to the underlying PLC+RDC contin- 



uum, we also detect a narrow emission feature consistent with neu- 
tral iron, and hence include a second, distant reflector. Here, as for 
NGC 7469, when allowing Cpin /xis to vary we were not able to 
reliably constrain all the reflection parameters of interest, so we 
present the results with the inclination fixed at 45° . 

Mrk 359: In addition to the basic PLC+RDC continuum, we 
also detect a narrow emission feature consistent with neutral iron, 
and hence include a second, distant reflector. We also include an un- 
resolved Gaussian component at ~5.9keV to account for a feature 
in the background spectrum. The formal global minimum gives a 
solution in z > 85° and the spin is unconstrained. However, as with 
PKS 0558-504 and IH 0323+342, we consider this to be unphysi- 
cally high for an unobscured source, and there is a local minimum 
of essentially equal statistical quality with i ~ 45° , so we instead 
present this solution. The spin is only weakly constrained. 

MCG-2-14-9: In addition to the PLC+RDC continuum, we 
also detect a narrow emission feature consistent with neutral iron, 
and hence include a second, distant reflector. A very good fit is 
obtained, however in this case, owing to the moderate data quality 
and the relatively weak reflection features, even when fixing the 
inclination at 45°, the iron abundance at the solar value, and the 
emissivity index at the Newtonian prediction (g = 3) we are unable 
to provide a preliminary spin constraint. 

ESO548-G081: In addition to the basic PLC+RDC contin- 
uum, we also detect a strong narrow emission feature consistent 
with neutral iron, and include a second, distant reflector. As with 
MCG-2-14-9 a very good fit is obtained, but owing to the moder- 
ate data quality and the relatively weak reflection features, even fix- 
ing the inclination at 45° and the emissivity index at the Newtonian 
prediction we are unable to provide a preliminary spin constraint. 

Mrk 1018: In this case, the PLC+RDC continuum provides an 
excellent fit, no statistically compelling narrow lines are detected. 
However, the spin is only very weakly constrained. 

RBS1124: Here, the basic PLC+RDC continuum also pro- 
vides a good fit to the available data . The s pin co nstraint for this 
dataset initially presented by Miniutt i et alj | |20IO|) is broadly con- 
sistent with that obtained here. 

IRAS 13224-3809: In this case, the XIS data quality is fairly 
poor, only extending up to ~7.6keV, and even relaxing our binning 
constraint to a S/N of 2 we do not obtain a robust PIN detection. 
However, we have recently undertaken a long observing campaign 
on this source with XMM-Newton, in which the source is found to 
display a sharp, deep drop above ~7.5 keV which is not well de- 
tected h ere owing to the poor data quality (s ee lFabian et alj|2012al . 
and also lPonti et alJl201()l . lBoller et alj 20031) . This feature, and the 
spectrum at higher energies, is extremely important in determining 
the iron abundance, inclination and also the photon index, so in this 
case we take these values from our XMM-Newton analysis. Further- 
more, to remain consistent with this more detailed analysis, we also 
include a low temperature (kT — 0.106 ± 0.005 keV) blackbody 
component, which may also be associated with irradiation of the 
accretion disc. In addition, we also include an unresolved Gaussian 
at ~5.9 keV to account for a feature in the background spectrum. 

Owing to the extreme nature of this source, by fixing the above 
parameters we still obtain a very strong spin constraint despite the 
relatively poor data quality, which is consisten t with the inner ra- 
dius originally obtained by IPonti et aL I ( l2010l) . and the spin ob - 
tained with our long XMM-Newton exposure jFabian et al.ll20I2^ . 
This remains the case even when the key parameters informed 
by our XMM-Newton analysis are allowed to vary within their 
quoted uncertainties. If we allow instead for a dual-ionisation re- 
flector (^1 — 20I7, = 490^200 ^•'g '-ni s~^), 3S proposed in 
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Figure 2. confidence contours for the spin measurements obtained for tlie compiled sample, based (where possible) on the results obtained when allowing 
Cxis/piN to vary (see Table[3). Horizontal dashed lines represent Ax^ equivalent to the 90 and 99 confidence intervals. 
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iFabian et al. I( l2012ij) . the fit improves to = 417/410, but the spin 
constraint obtained does not change. Attempts to fit lower spins, 
while still imposing the same limits on the other key parameters, 
result in the model significantly underpredicting the data at high 
energies. Removing the constraint on the photon index results in 
r — 2.22j;g Q5 and a substantially weaker constraint on the spin, 
a* > 0.9, although this is still constrained to be high. If the con- 
straints on the iron abundance and the disc inclination are also re- 
moved, the model becomes too degenerate to reliably constrain the 
spin solely using the Suzaku data, even if the blackbody compo- 
nent is also removed, and the values obtained for these quantities 
are substantially different to those obtained with the significantly 
higher quality XMM-Newton data (Ape/solar = 5^1, i = 30 ± 4°), 
despite not being expected to evolve with time. Furthermore, we 
stress that neither of these cases with relaxed parameter constraints 
correctly include the sharp high energy drop known to be present 
in this source. Therefore, we will proceed with the model in which 
the photon index, iron abundance and disc inclination are fixed at 
the values obtained wilh XMM-Newton . 

IH 0707-495: Again, the XIS data quality here is rather poor, 
the source only being detected up to ~6.7keV, and there is no 
robust PIN detection. As with IRAS 13224-3809, this source is 
well documented a s having a deep, sharp drop above ^7keV 



I Boiler et al 



IZoghbi etal 



.20021: iGallo et al 
l2010| - |Fabian et al 



20041: iFabian et alj |2004 |2009| ; 



l2012bh . which is not detected ow- 



ing to the poor data quality, leading to the same problems in de- 
termining the iron abundance, disc inclination and photon index. 
The spectrum obtained with this Suzaku observation is fairly rem- 
in iscent of the most re cent XMM-Newton observation, presented 
in IFabian etalj ( l2012bl) . so we adopt the photon index of 2.7 ob- 



tained in that work. We also fix the inclination and the iron abun- 
dance val ues at those obtaine d in the more detailed studies pre- 
sented bvlZoghbi etal.U2010l) | Fabian et al. and most re- 
cently l^aeTaiJIioii) . 



As with IRAS 13324-3809, by fixing these key parameters 
to their previously obtained values we are still able to obtain a 
strong spin constraint despite the poor data quality, owing again 
to the extreme nature of this source. The spin obtained here i s con- 
sistent with that originally presented in IFabian et al. I ( l2009l) . and 
also that found in subsequent reflection-based analyses (see also 
iDauser et alj|2012l) . Again, allowing the key parameters informed 
by the higher quality XMM-Newton data to vary within their quoted 
uncertainties does not significantly effect the spin constraint ob- 
tained, nor does adopting a dual-ionisation reflecto r (^i < 41, 
^2 = IS Otfin" erg cm s~^) similar to that proposed in lFabian et al.l 
( l2012d) . despite the latter modification providing a moderate im- 
provement in the quality of fit (xt = 264/234). Attempts to force 
lower spin values with the same constraints on the other key param- 
eters result in severe under-estimation of the data at higher ener- 
gies. In fact, this source was caught in such a reflection dominated 
state (see section l4~2l ) that, even removing the constraints on the 
photon index, iron abundance and disc inclination, the spin is still 
constrained to be high (a* > 0.88). However, the values obtained 
for the iron abundance and disc inclination (ylpe/solar = 3.41q g, 
i = 4015q°), which should not vary with time, are highly dis- 
crepant with those obtained from any high quality XMM-Newton 
observation, and as with IRAS 13224-3809 the sharp high energy 
drop known to be present disappears from the model. We therefore 
again proceed with the model in which the parameters are fixed 
based on previous XMM-Newton analysis. 



© 0000 RAS, MNRAS OOO.fTlllT] 



Suzaku Observations of 'Bare ' AGN 1 3 



lOl • r 



00 

o 5 - 



0.9 0.95 1 1.05 1.1 1.15 1.2 

X" / Degrees of Freedom 

Figure 3. The number distribution of reduced values obtained with our 
PLC+RDC interpretations for the presented sample of AGN. 

IRAS 05262+4432: Here, the XIS data quality is extremely 
poor, only extending up to ~7.8keV, and as with IRAS 13224- 
3809 and IH 0707-495 there is no robust PIN detection. A good 
fit is obtained with the basic PLC-fRDC interpretation, but even 
when fixing the inclination at 45°, the iron abundance at the solar 
value, and the emissivity index at the Newtonian prediction we are 
unable to provide a preliminary spin constraint. 



4 DISCUSSION 

We have presented a spectral analysis of a sample of 25 'bare' ac- 
tive galaxies observed with the Suzaku satellite. These are sources 
that display either no intrinsic absorption, or only very weak intrin- 
sic absorption, and have been selected on the basis of the presence 
of a relatively smooth soft excess below ~2keV after extrapola- 
tion of the ~2-10keV continuum to lower energies. Owing to this 
lack of absorption, such sources offer the best opportunity to study 
the intrinsic emission spectra of active galaxies. In addition to the 
soft excesses required for selection, a number of the sources also 
show evidence for broad iron emission components (e.g. Mrk335, 
Fairall9, Ark 564, see Fig. [TJ, and/or for 'hard' excesses above 
~IOkeV in the PIN data, although these features are not ubiqui- 
tous within the presented sample. 

Our analysis is based around interpreting these spectral com- 
plexities as reflected emission from the accretion disc. For each 
source we have constructed a continuum model consisting of a 
powerlaw-like component (PLC), associated with Compton up- 
scattering in an X-ray emitting corona, and a reflection dominated 
component (RDC) arising through irradiation of the surface of the 
accretion disc by the Comptonised emission. The reflected emis- 
sion has been modell ed with th e physically self consistent reflec- 
tion code REFLIONX jRoss & F abian 2005), which treats both the 
backscattered continuum emission and the atomic features pro- 
duced when optically thick material is irradiated by a powerlaw- 
like X-ray continuum, and the relativistic effects relevant for emis- 
sion from an accretion disc in the strong gravity regime close to 
a black hole are in cluded with the RELCONV convolution kernel 
jPauser et al]|2010l) . In addition to the their continuum emission, 
many of the AGN analysed also display narrow emission lines con- 
sistent with neutral iron, which we associate with reflection from 



cold, distant material, and treat with a second, unmodified RE- 
FLIONX component. A small minority also show evidence for dis- 
crete emission lines associated with highly ionised iron (Fe XXV 
and/or XXVI), which we model with simple Gaussian emission 
line profiles. All the emission components are modified by neutral. 
Galactic absorption. Finally, a further minority also show evidence 
for some minor absorption, either at low energies from moderately 
ionised material or at high energies from highly ionised material. 
Where present, we ha ve modelled this absorption with the XSTAR 
photoionisation code jKallman & Bautista 2001,) . In this work, we 
have attempted to simultaneously combine a systematic approach 
to analysing the sample as a whole, with a detailed physical treat- 
ment of each individual source (see section [T2] l. 

In all the cases analysed, we have been able to successfully 
model the available broadband spectra with a reflection based in- 
terpretation; Fig. |3] shows the distribution of the values obtained 
for the reduced chi-squared (xt = x^/D o.F.) with our PLC-l-RDC 
based interpretations for the analysed sample. This distribution 
clearly displays a fairly narrow peak around xt ~ ^- This demon- 
stration of the flexib ility of disc reflectio n is important, as we have 
previously argued in IWaltonetal.l ( l20I2h that this process is a fun- 
damental consequence of the widely accepted disc-corona accre- 
tion geometry. 

4.1 Sample Properties 

Having successfully modelled the observed data, we now consider 
the results obtained for some of the physical properties of the accre- 
tion flow for the compiled sample, focusing conservatively on the 
results obtained with Cxis/pin allowed to vary (Table [3). In gen- 
eral, when they could be constrained, the disc inclinations obtained 
from our analysis are either moderate or low, as would be expected 
for a sample selected on the basis of a lack of obscuratiorQ For a 
couple of the AGN considered the formal best fit does return un- 
physically high inclinations when this parameter is left free to vary 
given the nature of the selected sources. However, we attribute this 
to the complex dependencies and degeneracies that can arise be- 
tween the various parameters in such a complex, multi-parameter 
model, which can be especially problematic when modelling data 
of only moderate quality and/or a relatively featureless spectrum. 
Crucially, we are still able to obtain acceptable solutions in these 
cases even when fixing the inclination at a physically more reason- 
able value (45°). 

It may at first seem surprising that the majority of inclinations 
obtained are moderate (^-^50°) rather than low (<30°). However, 
there are two factors that result in a natural bias towards higher in- 
clinations (up to the point that sources are excluded due to obscu- 
ration by the torus, above ~70°). We have already mentioned that 
for disc-like structures, P{i) oc sin(i), so we naturally expect there 
to be more AGN with high inclinations. In addition, the majority of 
sources analysed appear to require high spin (see section |431 l. For 
such sources, the reflected flux observed in the soft X-ray band- 
pass is greater at high inclinations than low inclinations, owing to 
the larger line-of-sight velocities of the material in the disc at high 
inclinations. This results in an enhancement of the reflected flux 
observed owing to the additional special relativistic beaming along 
angles close to the plane of the disc, and also in additional reflected 
flux being blueshifted into the soft X-ray bandpass. The combined 

Recall that the probability distribution for the inclination of randomly 
orientated disc-like sources follows P{i) oc sin(i) 
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Figure 4. The number distribution obtained for the logarithm of the ioni- 
sation parameter of the reflecting medium with our PLC+RDC based inter- 
pretations. 
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Figure 5. The number distribution obtained for the logarithm of the iron 
abundance of the reflecting medium with our PLC+RDC based analysis. 



strength of these effects depends on both the spin of the black hole 
and the emissivity profile; for a maximally rotating black hole with 
a steep (g = 6) emissivity, the reflected flux at soft energies for a 
source observed at i = 80° can be more than a factor of ~2 larger 
than that for the same source observed at i = 10°, while for a non- 
rotating black hole with an emissivity of g = 3 the reflected flux in 
the soft band appears practically isotropic. Therefore, by basing our 
sample selection on the presence of a soft excess we are probably 
introducing an additional bias towards higher inclination sources, 
since rapid spins are generally inferred. 

In terms of the ionisation state of the reflecting medium 
{i.e. the surface of the disc), as shown in Fig. |4l we find a 
fairly broad range of ionisation parameters, ranging from ~10- 
1000 erg cm s~^. The vast majority of sources are found to have 
C < 400, suggesting that, as expected, iron is very rarely highly 
ionised in AGN. There is a fairly strong peak in the obtained dis- 
tribution around ^ ~ 300, co rresponding to iron b eing primarily 
found in its Fe XIX-XXI ions l lKallman et al]|2004l) . At these ion- 
isation states, resonant trapping and Auger destruction of the Ka 
photons is expected to be very efficient, and the emitted iro n line 
should only be very weak jMatt et alJl9"93l : lRoss et aljl996l) . This 
offers a natural explanation as to why a number of the sources do 
not show any strong evidence for relativistically broadened iron 
emission (see Fig.[TJ. 

In Fig. [5] we also plot the number distribution of the best fit 
iron abundances obtained; only sources for which we were able to 
constrain this quantity are included. The distribution clearly peaks 
around a solar iron abundance, with the vast majority of sources 
consistent with having an abundance within a factor of ~2 of this 
value. This is simi lar to the prefe r ence f or a roughly solar abun- 
dance reported by ICrummv et alj | |2006|) in their earlier applica- 
tion of disc reflection to a general sample of AGN, despite the 
two samples only sharing 5 sources in common. Our analysis re- 
sults in a substantially super-solar iron abundance for only one 
source, PDS 456 (although highly super-solar iron abundances are 
also adopted for IH 0707-495 and IRAS 13224-3809, this is based 
on prior analyses rather than being constrained by the data analysed 
here). As discussed previously, this is driven by the association of 



the emission feature observed at ~0.9 keV with iron L shell emis- 
sion in the reflected contribution. 



4.2 Reflection Strength 

A further quantity of interest is the implied strength of the reflected 
emission from the disc relative to the intrinsic PLC emission, i?disc, 
which can provide information on the geometry of the accretion 
flow. In order to estimate Rdisc, in a manner consistent with the lit- 
erature, we make brief use of t he PEXRIV Compton reflection code 
jMagd ziarz &Zdziarskilll995h . in which R is directly included as 
a free parameter. R is normalised such that a reflecting medium 
subtending a solid angle of SI = 27r, as seen by the illuminating 
source, gives a value of R= 1. The other key parameters of this 
model are the photon index and the high energy cut-off of the ion- 
ising continuum, and the inclination, abundances and ionisation of 
the reflecting medium. However, as PEXRIV does not self consis- 
tently include any flourescent atomic emission, which can be of 
significant importance in determining the reflected flux at soft en- 
ergies, we estimate i?disc from the strength of the Compton reflec- 
tion hump in our best fit models, rather than applying PEXRIV to 
the data directly, as follows. 

To remain consistent with our detailed modelling, the photon 
index, disc inclination, iron abundance and ionisation parameter 
are all set to the best fit values presented previously, and to re- 
main consistent with the manner in which the REFLIONX model 
was computed the high energy cut-off was set to 300 keV and the 
other elemental abundances were set to their solar values. The nor- 
malisation of the PEXRIV component was then set such that, when 
i?disc is set to zero the 2-10 keV flux matches that of the PLC com- 
ponent. Finally, 7?disc is determined such that the 15-50keV flux 
matches that of the combined PLC-l-RDC components in the rest 
frame, i.e. once the relativistic effects applied by RELCONV have 
been removed. We repeat this process for the best fit models ob- 
tained both when Cpin /xis was fixed and when it was free to vary, 
and the results obtained are given in Table|4l As we are not directly 
applying PEXRIV to the data, we do not estimate any statistical un- 
certainties on the values obtained for J?disc, and stress that these 
should only be considered as approximate estimates given the na- 
ture in which they are determined. 
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Table 4. Approximate estimates for tlie relative strength of the reflected 
emission from the accretion disc, based on the strength of the best fit Comp- 
ton reflection hump (see text). 
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Source 



Rd\ti 





CpiN/xis fixed 


CpiN/xiS free 


Mrk509 


0.7 


0.5 


3C382 


0.05 


0.1 


Mrk335 


0.8 


0.8 


Fairall 9 


0.4 


0.4 


IH 0419-577 


1.5 


1.3 


Ark 564 


1.5 


1.5 


Ark 120 


1.6 


1.5 


3C 390.3 


0.03 


0.03 


PKS 0558-504 


1.0 


1.0 


NGC 7469 


0.8 


0.2 


MrkllO 


2.6 


1.2 


Swift J0501. 9-3239 


2.1 


3.7 


Mrk841 


0.3 


0.3 


Ton SI 80 


1.4 


1.4 


PDS456 


1.2 


1.2 


IH 0323-1-342 


0.5 


0.5 


UGC6728 


1.5 


1.4 


Mrk359 


1.0 


1.0 


MCG-2-14-9 


0.3 


0.3 


ESO548-G081 


0.1 


0.1 


MrklOlS 


1.3 


1.2 


RBS 1124 


3.3 


2.9 


IRAS 13224-3809 


16 




IH 0707-495 


275 




IRAS 05262-H4432 


1.3 





In Fig.[6]we show the distribution of the values for iidisc ob- 
tained, where possible, when Cpin /xis was free to vary. This dis- 
tribution shows a clear peak around -Rdisc ^ 1 (the median value is 
~1.2), as broadly expected for a central isotropic emission source 
illuminating a geometrically thin accretion disc. However, there are 
also some specific individual results worth drawing attention to. 
First of all, the reflection fraction inferred for IH 0707-495 is ex- 
tremely large. This is consistent with our previous assertion that the 
Suzaku data is reminiscent of the XMM-Newton observation ob- 
tained recently in 2011, which caught IH 707-495 in an almos t 
completely reflection dominated state (see iFabian et alj l2012bl) . 
The reflection fraction obtained for IRAS 13224-3809 is also large, 
consistent w ith that estimated from our recent long XMM-Newton 
observation jFabian e t al. 2012a). At the other end of the scale, it is 
interesting to note that both of the radio-loud AGN included in the 
compiled sample display extremely weak reflection features. This 
would be expected if the illuminating X-ra y corona were associ- 
ated with the base of the jet, as suggested bv lMarkoff etal.U2005l) 
(see alsolSeloborodov 199^). It would be of great interest to further 
test this scenario by investigating whether systematically weak re- 
flection features are a general trait of radio loud AGN with a much 
larger sample of such sources. 



3 

o 
S3 



-2-10 1 

log iodise 

Figure 6. The number distribution obtained for log -Rdis 
our best fit PLC-hRDC continuum models. 



estimated from 



4.3 Black Hole Spin 

In addition to generally testing the robustness of the disc reflection 
interpretation, during the course of this analysis we have endeav- 
oured to use the inferred contribution from the reflected emission 
to place initial constraints on the black hole spin, although this has 
not always been possible, and it should be noted that the constraints 
presented are model dependent. Where the spin could be well con- 
strained, rapid prograde spins are generally obtained, although the 
black holes are certainly not always inferred to be maximally ro- 
tating. Even when the spin is only relatively loosely constrained, 
some kind of rotation is generally implied. There are no sources 
that are strongly constrained to be Schwarschild black holes, al- 
though there are a few sources that remain consistent with this in- 
terpretation. There are also no sources that are strongly constrained 
to have a retrograde spin, although again there are a few sources 
that remain consistent with this scenario. Therefore, if relativistic 
disc reflection is responsible for the intrinsic spectral complexities 
displayed by AGN, and as discussed previously there is very good 
evidence to suggest this is the case, the implication is that the ma- 
jority of AGN may be rotating fairly rapidly. Where the spin has 
been successfully measured via reflection for other sources not in- 
cluded in this work, the results obtained also seem to be consistent 
with t his picture (Brenneman & Reynolds 2006; Bre nneman et al.l 
I2OIII : lReisetal.lt2012: : Cowperthwaite & Revnolds.20121 : note that 
in the latter case, although the favoured spin implies a retrograde 
rotation with respect to the material in the disc, the black hole is 
still inferred to be rapidly rotating). For the majority of the sources 
analysed, where either a measurement of the spin or of the inner 
radius of the disc has been made previously, either using one of 
the same datasets analysed here or a different observation obtained 
at some other epoch, the spin obtained here is consistent with the 
previous work, although there are some notable exceptions high- 
lighted in section [T2I Where it has not been possible to reliably 
constrain the spin, this has largely been due to either a data quality 
issue, and/or a lack of strong reflection features. 

In this work we have carefully attempted to determine the cor- 
rect cross-normalisation constant to use between the combined XIS 
and PIN spectra based on the recommendations of the Suzaku HXD 
calibration team. To do so, we have first compared the combined 
XIS spectrum with the single XISO spectrum, and used this com- 
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Figure 7. A comparison of the spin constraints obtained with Cxis /pin 
fixed and with Cxis /pin fr^s to vary. Where limits are obtained, these are 
indicated with arrows. In many cases, allowing Cxis /PIN t° ^^ry does not 
significantly change the spin constraint obtained, although as discussed in 
section IX2l there are some notable exceptions. 



parison to calculate the normalisation constant Cxis/pin from the 
normalisation constants obtained by the Suzaku team between the 
XISO and PIN detectors using equation[T] In general we find a good 
agreement between the combined XIS and XISO spectra, as would 
be expected, although in rare cases the normalisations of these 
spectra can be discrepant by up to ~8 per cent (see UGC6728). 
We have also investigated what effect allowing this value to vary 
has on the results obtained, as it does not account for any systematic 
uncertainties on the level of the modelled instrument background. 
These are estimated to be -^3 per cent, and so we have allowed 
Cxis/pin to vary within an appropriate range based on this un- 
certainty. The results obtained allowing Cxis/pin to vary in this 
manner are quoted in Table[3] Focusing on the spin parameters, do- 
ing so does not appear to have any single systematic effect on the 
results obtained, either in terms of the preferred value or its statis- 
tical uncertainty. As demonstrated in Fig.|7] allowing Cxis/pin to 
vary has little or no effect on the spin obtained in many cases, even 
when the PIN spectrum carries some significant statistical weight 
(e.g. Mrk841), while in some the spin becomes more tightly con- 
strained (e.g. PKS 0558-504), and in others it becomes substantially 
more difficult to constrain {e.g. Mrk509). However, the general 
conclusions drawn above still seem to be robust. 

In Fig.|9]we plot a fractional distribution of the results shown 
in Fig.[2]for the black hole spins of the compiled sample. Here, we 
again only consider sources for which the data quality and/or the 
strength of the reflected features were sufficient to allow at least 
some constraint to be placed on the spin at the 90 per cent con- 
fidence level. However, in many cases the statistical uncertainties 
obtained are non-Gaussian and cover a broad range of spin val- 
ues, which can be significantly larger than any desirable bin size 
for such a distribution. As a compilation of AGN spin is one of 
the primary goals of this work, we take a more careful approach 
and attempt to account for these uncertainties when plotting the 
distribution of spins we obtain. We do so by combining the likeli- 
hood distribution^ for the spin parameters, L{a*), of the sample 




Figure 8. The likelihood distribution obtained for the spin parameter for 
Mrk335 (bottom panel) from the confidence contour {top panel); each 
bin spans a range of Aa* = 0.1. 



rather than simply the best fit values. These distributions are de- 
rived from the confidence contours displayed in Fig.|2l assum- 
ing that L(a*) cx exp(— and normalised such that the sum 
of the likelihoods over the —0.998 ^ a* ^ 0.998 range considered 
here is unity. As an example, in Fig.[8]we show the likelihood distri- 
bution obtained from the x contour for Mrk 335. Finally, the indi- 
vidual likelihood distributions are then combined and renormalised 
again to form a composite likelihood distribution for the compiled 
sample. The distribution obtained clearly displays a strong peak at 
higher spin value s. This is consi s tent w ith the typically small inner 
radii obtained bv lCrummv et al.l ( l2006l) . 

As discussed previously, the spin of the supermassive black 
holes powering AGN is widely expected to be determined by their 
growth history. Over the lifetime of these nuclear black holes, 
galaxy mergers and prolonged accretion can impart enough angular 
mome ntum to determine t he spin distribution of AGN in the current 
epoch. ISerti & Volonteril | |2008|) calculated the local spin distribu- 
tions expected for various different SMBH growth scenarios un- 
der the assumption that the SMBH seeds were not rotating. When 
SMBH growth is dominated solely by mergers, i.e. galaxy merg- 
ers resulting in the coalescence of the two central SMBHs, the ex- 
pected spins appear to be fairly randomly distributed, with a slight 
peak at low spin owing to their assumption regarding the initial 
seeds. However, prolonged accretion is found to be the dominant 
mechanism in determining the expected distribution of AGN spins. 
When the accretion is ordered, the central black holes are generally 
spun up, leading to a strong preference for rapidly rotating black 
holes. If instead the accretion is chaotic, i.e. the angular momenta 
of the accreted material and the AGN are not typically aligned, the 
AGN will experience the opposite effect and be spun down, leading 
to a strong preference for slowly rotating black holes. Therefore, 
measuring the spin for a large sample of AGN has important impli- 
cations regarding the nature in which their black holes grow. Our 
observational analysis suggests a preference for at least some rota- 
tion, and probably for rapid rotation, which if taken at face value is 



These distributions are designed to represent the underlying probability 



distributions for the spin parameters given the available data, but as this is 
not a fully Bayesian analysis they cannot formally be classed as such. 
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Figure 9. Fractional (left panel) and cumulative {right panel) distributions for the black hole spins obtained for the compiled sample of AGN (see text); each 
bin spans a range of Aa* = 0.1. The obtained distribution is clearly peaked at higher spin values. 



most consistent with the scenario in which SMBH growth is domi- 
nated by ordered accretion. 

However, before drawing any firm conclusions, we must first 
consider whether the general high spin preference obtained here is 
likely to be repre s entativ e of the true underlying spin distribution. 
iBrenneman et al. l|) suggest there may actually be an obser- 
vational bias towards AGN with high spin, based on the widely 
held expectation that the radiative efficiency of accretion onto a 
black hole increases with black hole spin, which should result in 
AGN with high spin being brighter on average than their low spin 
counterparts. Based on this, Reynolds et al. (2012) estimate that 
for a flux limited sample drawn from an AGN population with a 
flat intrinsic spin distribution, only ~30 per cent of the sources 
in the sample would have a spin of a* < 0.5, despite the true 
value being 50 per cent. Furthermore, as demonstrated in Appendix 
lAl we find that it is easier to statistically constrain higher spins, 
which is probably a natural consequence of the form of the re- 
lation between blac k hole spin and the radius of the ISCO (see 
iBardeen et al ][i97l). Given our method for compiling the sample 
likelihood distribution, this may also introduce some slight bias to- 
wards higher spins. Nevertheless, where we have been able to con- 
strain the spin well at the 90 per cent confidence level, rapid rota- 
tion {i.e. a* > 0.5) is almost always required, and this subset of the 
sample compiled is large enough that even given this probable ob- 
servational bias we might have expected to find at least one source 
that is constrained to have low spin if the intrinsic distribution does 
not truly peak at higher spins. However, we stress that the sample 
presented is certainly not a complete, flux limited sample, so the 
picture that appears to be growing could yet change. 

A further important consequence of our analysis, independent 
of the presence of any potential biases, regards the radio loud/radio 
quiet AGN dichotomy. It has previously been suggested that the 
observational differences between these classes of AGN are due to 
differences in the s pins of the black holes p owering these sources 
(Blandford 1990; Wilson & Colbert! 1 19951 ; iModerski et all 1 19981 ; 
[Garofalo et al.ll2O10l) . In this picture, radio loud AGN host black 
holes with high spin that can launch powerful jets via magnetic 
ext raction of the rotational ene rgy of the black hole, as proposed 
by IBlandford & Znaie while radio quiet AGN have low 

spin, and cannot launch jets via this method. However, our anal- 



ysis implies that, if the PLC+RDC interpretation adopted here is 
correct, then this cannot be the case, as the majority of AGN in- 
cluded in this work are radio quiet, yet are fou nd to have high spin. 
This i s not to say that the process outlined bv lBlandford & Znaiekl 
j 19771) is not present at all; when jets are present they may well 
carry more power for AGN with high spin than for AGN with low 
spin. However, black hole spin would not be sufficient by itself to 
solely determine jet production in AGN. This is similar to the ob- 
served behaviour of Galactic black hole binaries (BHBs), where the 
presence of jets appear s to be determined by the accretion state of 
the black hole (see e.g. iFender et al.ll2004l) ; the same BHB can be 
observed both with and without jets at different epochs, despite the 
fact that its spin cannot have evolved substantially. Black hole spin 
may yet prove to be important for certain jet properties, but it does 
not appear to play a substantial role in determining whether jets are 
produced in the first place. 

Although we have presented the largest compilation of AGN 
spins to date, in order to build a more robust picture the number 
of sources with reliable spin measurements must still be increased. 
There are two important aspects to addressing this issue. First and 
foremost, as we have stressed previously, in many of the cases in 
which we have been unable to provide an initial spin constraint, 
this has been due to the data currently available having a relatively 
poor S/N. This deficiency can naturally be addressed with further 
long observations with existing instrumentation. The other major 
uncertainty on some of the individual constraints obtained in this 
work is the cross-normalisation of the XIS and PIN instruments. 
In order to address this issue, it will be vital to obtain observations 
of nearby active gala xies with the upcoming NuStar observatory 
l lHarrison et al.l201oh . which is due to begin operation imminently. 
NiiStar will offer continuous spectral coverage over the 5-80 keV 
energy range, which includes a key region of overlap with the 
soft X-ray spectrometers currently available, and being an imag- 
ing spectrometer will allow simultaneous measurement of the back- 
ground emission. Simultaneous NuStar and siihsr XMM-Newton or 
Suzaku observations will be able to provide spin measurements ro- 
bust to the current cross-calibrational uncerta inties. In addition, the 
microcalorimeter due to fly aboard Astro-H l lKokubun et alJllOO^ ; 
2014) will enable the profiles of the soft excesses and the iron lines 
to be resolved in much greater detail, enabling us to further test 
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the robustness of the reflection interpretation and improve the spin 
measurements obtained. Finally, the effect of any high spin bias 
will ultimately need to be tested in the future with large, volume 
limited samples of AGN. 



5 CONCLUSIONS 

We have investigated the relativistic disc reflection interpretation 
for a large sample of 25 'bare' active galactic nuclei, sources with 
little or no complicating intrinsic absorption, observed with Suzaku. 
In constructing reflection based models for this sample, we have 
taken the unique approach of attempting to simultaneously under- 
take a systematic analysis of the whole sample as well as a detailed 
treatment of each individual source, and find that disc reflection 
has the required flexibility to successfully reproduce the broadband 
spectra of all of the sources considered. Where possible, we also 
use the reflected emission to place constraints on the black hole 
spin for this sample of sources. Our analysis suggests a general 
preference for rapidly rotating black holes, which if taken at face 
value is most consistent with the scenario in which SMBH growth 
is dominated by prolonged, ordered accretion. However, there may 
be strong observational biases towards AGN with high spin in the 
compiled sample, which limits our ability to draw strong conclu- 
sions for the overall population at this stage. Our analysis also im- 
plies that the radio loud/radio quiet AGN dichotomy is not related 
to black hole spin. 
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APPENDIX A: MULTI-PARAMETER REFLECTION 
MODELS 

One of the major aspects of this work is attempting to constrain 
the black hole spin from the full contribution of the reflected emis- 
sion from the accretion disc observed in AGN, i.e. from the com- 
bination of the soft excess, any broadened iron K emission and 
the high energy reflected continuum. In order to do so, we make 
use of a physically s elf consistent reflection model, REFLIONX 
( lRoss&Fabianll2005h . combined with the latest relativistic con- 
volution kernal, RELCONV l lDauser et alj|201ol) . This combination 
results in a complex, multi-parameter model for the reflected emis- 
sion observed, including the ionisation of the surface of the disc, 
its iron abundance, the the radial emissivity profile, the disc incli- 
nation and of course the black hole spin. The combination of the 
ionisation parameter and the iron abundance in particular can have 
a strong effect on the broadband spectral shape. In order to simul- 
taneously constrain all of these parameters, some balance between 
the strength of the reflection features and the quality of data avail- 
able is naturally required. In this appendix, we present a series of 
simulations demonstrating that reasonable quality Suzaku data can 
provide reliable constraints for all of these parameters. 

Using the FAKEIT command in XSPEC, we simulate Suzaku 
data based on our underlying PLC-fRDC continuum model (mod- 
ified by neutral absorption), using standard response and back- 
ground files for the (front illuminated) XIS and PIN detectors. 
These spectra are rebinned to the same S/N requirements as the 
real data presented in the main body of the paper, i.e. a minimum 
S/N per energy bin of 5 for the combined FI XIS spectra, and 3 
for the PIN spectra. The 2-lOkeV flux for the PLC component in 
these models is set to 2 x 10~^^ erg cm~^ s~^, and for simplic- 
ity both the PLC and RDC components are required to contribute 
the same flux over the much broader 0. 1-1000.0 keV energy range 
(i.e. the energy range over which the version of the REFLIONX 
model utilised here is calculated). The exposure time for the sim- 
ulations is set to 125 ks, such that we are simulating a fairly good 
Suzaku observation of a fairly bright AGN, relative to the sample 
compiled. We repeat this process for a variety of input parameter 
combinations, varying the ionistion of the disc, the iron abundance 
and the black hole spin. The other key parameters adopted are A'^h 
= 3 X 10^° atom cm~^, T = 2,q = 6,i = 45°, relatively typical of 
the results obtained in the main body of the paper. For illustration, 
Fig. lAll shows some examples of the ratio plots obtained with the 
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Energy (keV) Energy (keV) Energy (keV) 

Figure Al. Data/model ratio plots for a selection of the simulated Suzaku XIS (front illuminated; black) and PIN (red) spectra to a simple powerlaw continuum 
model. The plots for the two higher spin values shown display reasonable qualiatative similarities with many of the real datasets shown in Fig.[T] The simulated 
data shown have been re-binned for plotting purposes only. 



Table Al. The various input parameter combinations simulated, and the values obtained modelling these simulated spectra with our PLC+RDC continuum; 
all the spectra are simulated with F = 2.0, g = 6 and i = 45° . The input parameters are generally successfully reproduced. 
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simulated spectra when modelled with a simple absorbed powerlaw 
model, applied as described in section [2^ The plots obtained for 
the two higher spin values, although obviously showing idealised 
scenarios in which just the PLC and disc reflection components are 
present, are reasonably similar to many of those obtained with the 
real data (see Fig.[Tll. Each of the simulated spectra are then mod- 
elled with our PLC+RDC continuum model, and the parameters 
obtained are presented in Table I All 



In general, the input parameters are successfully reproduced, 
particularly the parameters that significantly help determine the 
broadband spectral shape, i.e. the iron abundance, disc ionisa- 
tion and photon index. The relativistic blurring parameters are also 
generally well reproduced, although there is one interesting trend 
worth briefly discussing. It is clear from Table lAl] that it is easier 
to constrain the spin when the spin is high. This can be understood 
in terms of the evolution of the theoretical relation between the ra- 
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dius of the ISCO with black hole spin (see lSardeen et al J 19721) . At 
higher spins, a given deviation in spin results in larger deviations 
in the radius of the ISCO. As the model utilised here ultimitely in- 
fers the spin by determining the inner radius of the disc from the 
observed spectrum, and assumes this is coincident with the ISCO, 
it is therefore naturally easier to statistically constrain higher spins 
via this method. Nevertheless, these simulations demonstrate that 
Suzakii data with reasonably good S/N (as well as reasonably strong 
reflection features) should be sufficient to determine even moderate 
spins, as well as the other key parameters of interest in this work, 
from the observed broadband X-ray spectrum. 



Suzaku Observations of 'Bare ' AGN 
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